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Abstract 



^4 In this paper we obtain the density function and the distribution 

function of the distance between two uniformly and independently 



distributed random points in any right-angled triangle. The density 
function is derived from the chord length distribution function using 
Piefke's formula. We conclude results for random distances between 
two congruent right-angled triangles together forming a rectangle. 
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1 Introduction 

A random line g intersecting a convex set fC in the plane produces a chord 
of /C. The length s of this chord is a random variable. The chord length 
distribution function of a regular triangle was calculated by Sulanke |141 
p. 57]. Harutyunyan and Ohanyan [8] calculated the chord length distribu- 
tion function for regular polygons (see also Basel [2j pp. 2-9]). 

The distance t between two points chosen independently and uniformly 
at random from K, is also a random variable. Borel [3] considered random 
distances in elementary geometric figures such as triangles, squares and so 
on (see [HI p. 163]). The expectations of the distance between two random 
points for an equilateral triangle and a rectangle are to be found in [TBI 
p. 49]. Ghosh [6] derived the distance density function for a rectangle. Basel 
[2J pp. 10-16] derived the density function and the distribution function of 
the distance for regular polygons from its chord length distribution function 
using Piefke's formula |121 p. 130]. There are a lot of results concerning the 
distance t within a convex set or in two convex sets, see Chapter 2 in |10j . 

For practical applications in physics, material sciences, operations re- 
search, communication networks and other fields see [7], [9] and [15j . 

The first aim of the present paper is to derive the chord length distribu- 
tion function for any right-angled triangle. The second aim is to conclude 
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the distribution function of the distance between two points chosen inde- 
pendently and uniformly from a right-angled triangle. 

A = (0, b), B = (a, 0) and C = (0, 0) are the vertices of the right-angled 
triangle (see Fig.[[b). We assume \AC\ = b > a = \BC\ and denote this 
triangle by T a ,b- Clearly, it is possible to bring every right-angled triangle 
in such a position by use of Euclidean motions (translations and rotations) 
and, if necessary, a reflection. These operations do not influence the chord 
length distribution and the point distance distribution, c = \AB\ is the 
length of the hypotenuse and a = arctan(a/6) the inner angle between AB 
and AC. u = a + b + c is the length of the perimeter of T a ,bi A = ab/2 its 
area and h = a cos a the height over AB. The hypotenuse AB of the triangle 
Ta,b is the diameter of its circumcircle (converse of Thales' theorem). 



2 Chord length distribution function 

A straight line g in the plane is determined by the angle ip, < ip < 2ir, 
that the direction perpendicular to g makes with a fixed direction (e.g. the 
x-axis) and by its distance p, < p < oo, from the origin O: 

9 = QiPiip) = {{ x -,y) ^ : x cos ip + ysintp = p} . 

The measure /J, of a set of lines g(p, ip) is defined by the integral, over the set, 
of the differential form dg = dp dip. Up to a constant factor, this measure 
is the only one that is invariant under motions in the Euclidean plane [13, 
p. 28]. 

We use the chord length distribution function of T a ,b in the form 

F{s) = - fj,({g : 9 nr a ,^0, \ X (g)\ < a}) , 
u 

where x(<?) = 9^Ta t b is the chord of 7^,6, produced by the line g, and |x(<?)| 
the length of x(g) (see [5, p. 863] and [1, p. 161]). (The measure of all lines 
g that intersect a convex set is equal to its perimeter [THJ p. 30].) So it 
remains to calculate the measure of all lines that produce a chord of length 
|x($)l ^ s ■ Using the abbreviation 

S(g,s) := {g: 9 nT a ,^0, \ X (g)\ < s} , 

we have 

K s (9, s)) = dg= dp dip. 

JS(g,s) JS(g,s) 

Instead of the angle tp it is possible to use the angle cp that a line makes 
with the positive x-axis, hence 



K s (9, a)) = 




dp d(p . 
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Theorem 1. The chord length distribution function F of T a ,b is given by 



F(s) 



( if s<0, 

Hi(s)/u if 0<s<h, 

H2(s)/u if h<s<a, 

Hs(s)/u if a<s<b, 

Hi(s)/u if b<s<c, 

1 if s > c, 



where 
Hi(s) 
H 2 (s) 



9i s 

e 2 s 



+ 6 3 Li(s, h) + cL 2 (s, h) , 



H 3 (s) = a + 6 4 s + Li(s,/i) + ^L 2 (s,h) + ^Li(s,a) + ^L 2 (s,a), 
H 4 {s) = 6 5 + ^ + ^ Li(a, a) + ^ L 2 (s, a) + |- Li(s, 6) + | L 2 (s, b) 



with 



2b 



0i = ^ (2a + tt) + — (tt - a) + 6 , 9 2 = £ (2a - vr) - — a + 6 , 
o a b a 

©3 = r + -> ©4 = 1 - - a , 6 5 = a + 6, 6 6 = ^ (2a - vr) - — a + 2 
o a a o a 



and 



a) = sarcsin - , L 2 (s, a) = \J 1 — ^—J . 



Proof. In the following we need the width w olT a ,b i n the direction perpen- 
dicular to the angle 4> that is given by: 

' W!(<f>) if < 4> < vr/2, 
w(4>) = I w 2 {<j)) if tt/2 < 4> < ir/2 + a , 
Wz{4>) if 7r/2 + a < <p < TT , 



where 



/ ,s cos(<p - a) . cos 

wi((f>) = a : , W2\<p) = asmq) , w^{(p) = —a 



sin a 



tan a 



We have to distinguish the cases < s < h, h < s < a, a < s < b and 
b < s < c, and put 

ffi(s) = n(S(g, s)) if 0<s<h, H 2 (s) = n(S(g, a)) if h<s<a, 
Hs(s) = n(S(g, s)) if a < s < b , H^s) = n(S(g, s)) if b < s < c . 
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We demonstrate the calculation of H 2 . Here it is necessary to distinguish 
the following intervals of (ft (see Fig. [T|) : 

a) O<0<0i=a — arccos ((a/s) cos a) : All lines g with fixed angle 
produce a chord x °f length |x(<50l < s if they are lying in one of the two 
strips of breadth 

, . . s sin(20 — a) + sin a , . s cos(20 — a) + cos a 

0=o and &2(s, 0) = o = 

2 cos a 2 sin a 

respectively. 

b) <ft\ < <ft < 4>2 = a + arccos ((a/s) cos a) : All lines g with fixed angle (ft 
and g n 7^,6 7^ produce a chord of length \x(g)\ < s. These lines are 
lying in a strip of breadth w\(4>). 

c) ^2 < (ft < vr/2 : See < (ft < <ft x . 

d) 7r/2 < < 7r/2 + a : All lines g with fixed angle produce a chord of 
length < s if they are lying in one of the two strips of breadth 6i(s, 0) 
and 



h(s,(ft) = -- sin 20 



respectively. 



e) 7r/2 + a < (ft < ir : All lines g with fixed angle 0, lying in one of the 
two strips of breadth 62(5,0) and 63(5,0) repectively, produce a chord 
of length < s. 

Therefore, using the abbreviations b{ and Wi for 6j(s,0) and u>i(0) respec- 
tively, 

-7r/2 



H 2 (s)= [ dpd(f>= [ 1 (b 1 + b 2 )d(/)+ [ * Wl d(ft+ [ (61 + 62) 

JS(g,s) JO J<f>i J(f>2 

+ / (6! + 6 3 )d0 + / (6 2 + 6 3 )d0. 

Jir/2 Jw/2+a 



With the angles 

, a . fa 

03 = 7T — arccos — and 04 = it — arcsm I — cot a 
s \s 

we analogously get 

Hi(s)= r /2 (6i + 6 2 )d0+ r /2+Q (6i + 6 3 )d0+ r (6 2 + 6 3 ) 

JO J-k/2 Jn/2+a 

r4>2 /"""/ 2 r-n/2+a 

H 3 (s) = Wl d4>+ (6 1 + 6 2 )d0+/ (6i+6 3 )d0 

JO J<f> 2 Jk/2 



/ (6 2 + 6 3 )d0+ / w^dcft, 

Jtr/2+a J<f>z 
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HA s) 



«r/2 






/ vox d 




<f> + / 


/o 


Jir/2 





n/2+a 



(h + 63) • 



r<t>3 fir 

+ / (b 2 + b 3 )d(j)+ / w 3 d0. 

Jir/2+a J 6* 



r/2+a J <f) 3 

Using the indefinite integrals 

s 2<f) sin a — cos(20 — a) 



J &iM) d ^ 

b 2 (s,(f>) dcp 



4 cos a 

s 2</> cos a + sin(20 — a) 

4 sin a 



63(5, 4>) d<j) = — cos 2 <j) . 



w 1 (4>) < 
J w 2 (<f))< 



asm [0 — a 



sin a 



= — acos^>, / ws(4>) deft 
one finds the result of Theorem [TJ 



a (cot a sin (ft — cos< 



a sin 
tana 



□ 



Remark. It is easy to derive the chord length density function from Theo- 
rem [TJ 

3 Point distances 

Theorem 2. The density function g of the distance t between two random 
points in T a b is given by 



9{t) 



2t 
1 




if 0<t<c, 
else , 



where 



r(t) 



( J*(o,t) 
Jt(p,h) + J$(h,t) 

Jf(0, h) + J$(h, a) + J 3 *(a, t) 



if 0<t<h, 
if h < t < a , 
if a <t < b , 



with 



k Jt(0,h) + J%(h,a) + J£(a,b) + JZ(b,t) if b<t<c 
J* k (s,t)=H* k (t)-H* k (s), 
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e^ 2 



where 
Hl(t) = at + 



- + Q 3 Ll(t,h) + cL* 2 (t,h), 

e 4 t 2 , e 3 



+ -f Ll(t, h) + - L* 2 (t, h) + — L\{t, a) + - L*(t, a) , 



with 



. a 
a arcsm — 
t 



L\(t } a) = - (a\Jt 2 - a 2 + t 2 arcsin |) , L%(t, a) = \Jt 2 - a 2 + 

Proof. According to a theorem of Piefke |12} p. 130], the density function of 
the distance between two random points is given by 



. . 2ut r 

alt)= pL 



«)/(«) d., 



where / is the density function of the chord length. Following [2, pp. 11/12], 
we obtain 



2t 



7T ~ [t 

A 



7T + — {J* {t) - Ut) 



with 



This yields for < t < h 
J*(t) 

and for h < t < a 



F{s)ds 

J*(t) := u I F{s)ds. 
Jo 

[ H 1 (s)ds = H{{t) - H*(0) = Ji(0,t) , 
J o 



J*(t) = u / F(s)ds 



[ H 1 (s)ds + [ H 2 (s)ds 

JO Jh 



= H{(h) - fl?(0) + fl|(t) - H* 2 {h) = Jf(0, /») + J 2 *0, i) . 

The calculations of J*(t) for the intervals a < t < b and b < t < c can be 
carried out in the same manner. 

Now we determine the indefinite integrals of L\(t, a) and L 2 (t, a). Using 
integration by parts, we find 



Ll(t,a) 



Li(t, a) dt 



— ( I" arcsin — + a 



t arcsin — dt 
t 



VW^ 2 



dt 
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Clearly, 



and therefore 



Vi 2 - a 2 



dt = Vt 2 - a 2 



Ll(t,a 

Since t > a > in the present cases 



— (ay/ 1 2 — a 2 + t 2 arcsin — 



L* 2 (t,a) 



L2 (t, a) dt 



dt 



dt. 



and (cf. H p. 48, Eq. 217]) 



L*(M) = V* 2 



a 1 + a arcsin • 



□ 



Corollary 1. The distribution function G of the distance t between two 
random points in l~ a ,b is given by 



G(t) 




1 

I 
1 



with 



J\t) 



r 

K x (h) + K 2 (a) + K 3 (t) 



if t<0, 
if 0<t<c, 
if t>c 



if 0<t<h, 
if h < t < a , 
if a < t < b , 



k K x {h) + K 2 (a) + K 3 {b) + K 4 (t) if b<t<c, 

where 

K^t) = 4(0,t) , 

K 2 (t) = \( t2 - h ") i J i (0. >0 " #2 CO] + jJ(M) , 

K 3 (t) = \{t 2 - a 2 ) [J*(0,/i) + J*(h,a) - H*(a)} + j\{a,t) , 



K, 



(t) = \{t 2 - b 2 ) [Ji(0,h) + JZ(h,a) + J 3 *(a,6) - Ht{b)\ + j\(b,t) 



with Jt and Ht according to Theorem and 



jl(s,t) = Hl(t)-Hl(s) 



where 
H\(t) 

H\{t) 



01 t A 

32 
9 2 t 4 

32 



+ 3 L}(t,/i) + c4(t,/i), 



H\{t) = ^ + ^ + ^ L \(t,h) + ^4(t,h) + ± L\(t, a) + \ 4(t, a) , 
H\(t) = ^f + ^f + l4(t,a)+ b - L\(t, a) + ± L \(t, b) + \ L>(t, b) 



with 



L\(t,a) 



^ (f _ a 2 f 2 + a 3 y^T^ + t 4 arcgin £ 



L 2 (t, a) = \ {t 2 - a 2 ) 3/2 + £ (ay 7 ^^ 2 + ^ arcsin 

3 2 \ i / 



Proof. For < £ < c one gets 



G(t)= f g(r)dr = f 
Jo Jo 



2txt 2ut 2 2ut 



A 



A 2 A 2 



[ F(s)di 
Jo 



s I dr 



TTt 2 


2ut 3 


2u 

+ u 


A 


3A 2 


TTt 2 


2ut 3 


2 


A 


3A 2 



1 

A 



3A 



A 



with 



j\t) := f sJ*(s)ds. 
Jo 

It remains to calculate J\t). For < t < h we have 

j\t) = I sJ${0, s)ds= f s [HZ(s) - HZ(0)] ds 
Jo Jo 

= H\(t) - tf}(0) - \t 2 H{(V) = H\(t) - H\(0) - = 4(0,t) = K^t) . 
For h < t < a we find 

j\t) = K 1 (h)+ [ s[J?{0,h) + J 2 *(/i,s)]ds 
Jh 

1 /"* 

= K 1 (/i) + -(t 2 -/ l 2 ) jr(o,/i)+ / sj;{h, s )ds. 

2 Jh 

With 

/ sJ^(h,s)ds = f s[H* 2 {s) - H* 2 {h)]ds = f sH^{s)ds - H%(h) f sds 

J h J h J h J h 
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H\{t) - E*(h) - I [t 2 - h 2 ) HI (h) 
j\(h,t)-\{t 2 -h 2 )Hl{h) 



it follows that 



J\t) = K x (h) + \{t 2 ~ h 2 ) [JUO, h) - H* 2 {h)\ + j\{h,t) = K^h) + K 2 (t) . 

The calculation applies analogously to the intervals a <t < b and b < t < c. 

Now we calculate the indefinite integrals of L*(t, a) and a )- We 
have 

L\ (t,a) = J tL*(t,a)dt = ^ J t (ay/t 2 - a 2 + t 2 arcsin ^ dt 
= ~(ajt \Jt 2 - a 2 dt + J t 3 arcsin j dt\ 

with 

J t\/t 2 -a 2 dt = i (t 2 - a 2 f' 2 @1 p. 47, Eq. 214] . 
Using integration by parts, we find 

f 3 ■ a , 1 ( 4 ■« /" * 2 , \ 

/ t arcsin — dt = - \ t arcsin — ha/ — = dt . 

y t 4 v * y vi - (a/*) 2 / 

Since t > a > in the present cases, 

t 2 .ft 3 



y/l - (a/t) 2 J y/t 2 - a 2 



dt = / , dt . 



and 

f3 



J J r _^ dt=\ (t 2 - a 2 f 2 + a 2 V¥^7 2 ® p. 48, Eq. 223] . 
This yields 



L\(t,a) = l 



^ [t 2 _ a 2 f 2 + a 3 ^T—I + t 4 arcgin « 
3 i 



Further, 



L^t, a)= y tL2(t,a)dt = J ty/ 1 2 — a 2 dt + a ^ t arcsin ^ dt 



= - (t 2 - a 2 ) 3/2 + ~ (a\/t 2 - a 2 + t 2 arcsin - 
o 2 V t > 

(see the calculations of l\{t,a) and L*(t, a)). □ 
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4 Two triangles 

Now considering a rectangle 1Z a ^ of side lengths a and b (a < b) as made 
up of two congruent copies of the triangle T a , b with different orientation and 
sharing their hypotenuses. Let g 2 and G 2 denote respectively the density 
function and the distribution function of the distance between two uniformly 
and independently distributed points Pi and P2, one in each of the triangles. 

Using the indirect method (see Ghosh [6j pp. 22/23]), it is possible to 
calculate g 2 from the density function gn of the distance between two ran- 
dom points inside TZ a ,b and the density function g of T a ,b (according to 
Theorem [2]). The probability of Pi and P2 belonging to the same triangle is 
equal to the probability of Pi and P 2 belonging to different triangles, each 
probability being 1/2. Therefore, 
1 



giz 



{g + 92) , whence g 2 = 2 g n - g and G 2 = 2 G n - G . 



where Gn denotes the corresponding distribution function of gn- As already 
mentioned, gn was found by Ghosh [6] (Theorem 2, pp. 18/19). Denoting 
by A-ji = ab the area of TZ a ,b and by u-ji = 2 (a + b) its perimeter, we can 
write Ghosh's result as 



gnit) 



A 2 ^(*) 





if < t < c 



v 7 ^ 



+ 6 2 , 



else , 



where 



H ni {t) 
H h..2( t ) 



l 1l,2\ 



if 
if 
if 



< t < a. 
a <t <b. 
b<t< c 



with 



ttAti - unt + t , 
-a 2 -2bt + 2bL* 2 {t } a) 

„2\ .2 



H n ^{t) = -(7rA n + c z )-t z + 2bL* 2 (t,a) + 2aL* 2 (t,b) 

(L 2 as in our Theorem [2]). From this, one easily concludes the distribution 
function 



G n {t) 



gn(r) dr 




1 



if 
if 
if 



where 



Jl 



A(t) 



J ni (0,a) + J^ 2 (a,t) 
4 1 (0,a) + 4 2 (a,6) + 4 i3 ( 



t<Q, 
0<t<c, 

t > c, 

if < t < a , 
if a <t <b, 
if b<t<c 
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with 
where 

H n,i^) = J tH n,i(t)dt = — — + - , 

/„2j.2 9hf3 
tH^ 2 (t)dt = -- - + 2bL\(t,a), 

< 3 (t) = f tH^ 3 (t) dt = - M * ± ^ - ~ + 26ij(t, a) + 2a4(t, b) 

with as in Corollary [TJ Now it is easy to calculate the density func- 
tion 52 and the distribution function G<i- Examples are to be found in 
subsection 15.21 



5 Examples and simulation 

5.1 One triangle 

Fig. [2] shows some density functions eg. The representation in this diagram 
has the advantage that eg is independent from the maximum distance c, 
it depends only on the ratio a/b. Fig. shows the comparison between 
distribution functions G and corresponding simulation results (plot markers) 
with c = 1 and 10 7 pairs of random points. 



C xg(t) 




0.2 0.4 0.6 0.8 1.0 ' 

Figure 2: T a , a (thick), T a ,2a (thick and dashed), T a ,5a (thin and dashed) and 
T a< 20a (thin) 
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Git) 




Figure 3: T a , a (thick) and T a ,5a (thin and dashed) 
5.2 Two triangles 

Fig. H] shows the density function cgi for b = 5a, and the corresponding 
functions cg-ji and eg. Fig. [5] shows the distribution function G2 for b = 5a, 
and the corresponding functions G-ji and G. The plot markers result from 
a simulation with c = 1 and 10 pairs of random points. 

c xg(t) 



3.0 - 




Figure 4: cgn (short-dashed), eg (long-dashed) and C02 for b = 5a 
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G(t) 



1.0 




0.8 



0.6 



0.4 



0.2 



t/c 



0.2 



0.4 



0.6 



0.8 



1.0 



Figure 5: G-ji (short-dashed), G (long-dashed) and G2 for b = 5a 
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